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Abstract 

We study the homogenization problem of semi linear reflected partial differential equa- 
tions (reflected PDEs for short) with nonlinear Neumann conditions. The non-linear term 
is a function of the solution but not of its gradient. The proof are fully probabilistic and 
uses weak convergence of associated reflected generalized backward differential stochastic 
equations (reflected GBSDEs in short). 
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1 Introduction 

Let L% be an uniformly elliptic second-order partial differential operator indexed by a parameter 
e > 0. The homogenization problem consists in analyzing the behavior as e j 0, of the 
solution u £ for the partial differential equations (PDEs): L £ x u £ — f in a domain of IR d . 
The coefficients of the operator L e x are subject to an appropriate conditions. For example, 
periodicity by Freidlin [5], almost periodicity by Krylov [9]. In this paper we shall consider 
the case of locally periodic coefficients introduced by Ouknine et al. [33]. In the probabilistic 
approach the problem becomes the following. What is the limit of the laws of the diffusion 
processes X s with generator L x , as e { ? This problem has been studied for diffusion processes 
in the whole of IR d by Krylov|9], Freidlin[5] and Bensoussanp]. In the case of the presence of 
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boundary conditions, we can see the work of Tanaka [18] and recently this one of Ouknine et 
al.[13]. In [15], Pardoux has combined the probabilistic approach of chapter 3 of [T] to derive 
homogenization results for semi linear parabolic partial differential equation involving a second 
order differential operator of parabolic type in the whole space with highly oscillating drift 
and nonlinear term. Furthermore Pardoux and Ouknine [12] pursue the same program for a 
semi-linear elliptic PDE with nonlinear Neumann boundary conditions. 

This paper is devoted to prove a similar result for obstacle problem of semi-linear elliptic 
partial differential equations (PDEs in short) with a nonlinear Neumann boundary conditions g 
satisfying moreover linear growth conditions. We use for this purpose the link between backward 
stochastic differential equations (BSDEs in short) and PDEs. 

The concrete formulation of BSDEs in the nonlinear case was first introduced by Pardoux 
and Peng [14] who proved existence and uniqueness of adapted solutions for these equations, 
under suitable conditions on the coefficient and the terminal value. They provide probabilistic 
formulas for solutions of semi linear PDEs which generalized the well-known Feynman-Kac 
formula. The interest of this kind of stochastic equations has increased steadily, since it has 
been widely recognized that they provide a useful framework for the formulation of many 
problems in mathematical finance (see [§]), in stochastic control and differential games (see [7]). 
A class of BSDEs called generalized BSDEs has been introduced by Pardoux and Zhang in [16j. 
This kind of BSDEs has an additional integral with respect to a continuous increasing process 
and is used to provide probabilistic formulas for solutions of semi linear PDEs with nonlinear 
Neumann boundary conditions. Further, in order to provide probabilistic formula for solutions 
of obstacle problem for PDEs with nonlinear Neumann boundary conditions, Ren and Xia [TZ] 
considered reflected generalized BSDEs. Our goal in this paper is to give a homogenization 
result of obstacle problem for semi-linear PDEs with nonlinear Neumann boundary conditions 
via a weak convergence of reflected generalized BSDEs. 

To describe our result more precisely, we first recall some classical notations that will be 
used in the sequel of this paper. Let C^(IR ) be the space of real continuously differentiable 
functions such that the derivatives are bounded and 9 be a regular convex and bounded domain 
in IR d . We introduced a function p G C b 1 (IR d ) such that p = in 6, p > in IR d \9 and 
p(x) = (d(x,Q)) 2 in the neighborhood of 9. Let us note that, since the domain 9 is smooth, 
it is possible to consider an extension ip G C 2 (lR d ) of the function d(x,d@) defined on the 
restriction to 9 of the neighborhood of <99 such that 9 and <99 are characterized by 

9 = {x G IR d : ip(x) > 0} and <99 = {x G IR d : ^(x) = 0} 

and for all x G dO, Vip{x) coincide with the unit normal pointing toward the interior of 9 
(see for example [ID], Remark 3.1). In particular we may and do choose p and i/j such that 
(Vij)(x),5(x)) < 0, for all x G IR d , where S(x) = Vp(x) which we call the penalization term. 
Now let us consider the following reflected semi-linear partial differential equations with Neu- 
mann boundary condition for e > 0): 
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min (u e (t, x) — h(t, x), 

(t, x) - L £ x u £ - f(x, u £ (t, x))) = 0, if (t, x) e IR+ x 0, 

< 

T £ x u £ (t, x) + g(x, u £ (t, x)) = 0, if (t, x) G IR + x <90, 

u £ (0,x) = l(x), x e 0, 

where is defined in fl 2 . 2 j) . Under suitable linear growth, monotonic and Lipschitz condi- 
tions on the coefficients /, g and I, we will show that u £ (t,x) the viscosity solution of the 
previous reflected PDEs converge, as e goes to 0, to a function u(t,x), the viscosity solution to 
the following reflected PDEs with nonlinear Neumann boundary conditions: 

min (u(t, x) — h(t, x), 

-f (t, x) - Liu - f(x, u(t, x))) = 0, if (t, x) e H + x 0, 

< 

r° x u(t, x) + g(x, u(t, x)) = 0, if (t, x) G IR+ x 90, 

u(0, x) = l(x), x E 0, 
where operators L x and T x are in the form (12.41) . 

The rest of this paper is organized as follows. In section 2 we recall without proof, some 
results due to Ouknine et al.[13] on the homogenization of the reflected diffusion in the bounded 
domain with locally periodic coefficients and some basic notations on generalized BSDEs. Sec- 
tion 3 is devoted to the statement and proof of our main result. Finally, in Section 4, we give 
a homogenization result for some PDEs with nonlinear Neumann boundary conditions. 



2 Locally periodic homogenization of reflected diffusion 

In this section, we state some results which were proved in Ouknine et al. [13], and which we 
will use in the rest of the paper. Let us recall that is a regular convex and bounded domain 
in H d ; (d > 1) defined as previous and let 

d 1 ^ ^ 

L l = ^2 a ij( x , ~) d A + ~^2 b i(x, -)di +^2ci(x, -)di (2.1) 

i,j i=l i=l 

and 

d 

K = X>Vtof)3i (2-2) 
i=i 
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be given (with e > 0), where di = d/dxi. We also set 



i=i 



and require the following: 



(H.l) L XjJ , is uniformly elliptic and the matrix a(x) — [a^j] is factored as a(x,y)a*(x,y)/2 

(H.2) The functions a : IR d x -> IR d , b : IR d x 6 -> IR d , c : IR d x 6 -> IR d are locally 
periodic (i.e, periodic with respect to second variable; of period 1 in each direction in 0). 

(i) Global Lipschitz condition: there exists a constant C such that for any ( = a, b, c, 
\\((x,y)-ax>,y')\\ <C(\\x-x>\\ + \\y-y>\\), Vx,x'elR d , M 'e6, 

(ii) The partial derivatives d x ((x,y) as well as the mixed derivatives d% y £(x,y) 
exist and are continuous, Wx G IR d , y G for ( = a, 6, c 

(m) The coefficients are bounded, that is, there exists a constant C such that for any 
{ ( = a,b,c \\C(x,y)\\<C, xeM d ,ye&. 

The differential operator L £ x inside together with the boundary condition Y x u = on <90 
determine a unique diffusion process X s in which we call the (L £ x , r^)-diffusion. Given a d- 
dimensional Brownian motion (W t : t > 0) defined on a complete probability space (Q, J 7 , IP), 
let (X e , G e ) be the unique solution with values in x IR + of the following reflected SDE: 

{ dXf = \b{Xl &)dt + c{Xl &)dt + a(Xf, ^)dW t + V^XDdGl t > 



XI G 0, G e is continuous and increasing such that, 
jJV^(X«)l Wee }dG; = 0, t>0 



(2.3) 



By requirement there exists a L x ^-diffusion on IR d with generator L x y and by y-periodicity 
assumption on the coefficients this process induces diffusion process U x on the <i-dimensional 
torus T d , moreover this diffusion process is ergodic. We denote by m(x, .) its unique invariant 
measure. In order for the process with generator L e x to have a limit in law as e [ 0, we need 
the following condition to be in force. 



(H.3) Centering condition: for all x, 



L 



b(x, u)m(x, du) = 0. 
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Let us denote 



\b d J 

with b k (x,u) = / TE u {b k (x, U*)}dt, where under JP U , U x starts from u. 
Jo 

We put 

d d d 

/-'; = E ^(^^ + E^)^ r° = E-:v)^, 

i,j=l i=l i=l 

where the coefficients are respectively defined by 

1 



<9 X 66 + (J + S^6)c + -Trdiybaa 



(2.4) 



Sq{x,v)= (I + d y b)a , A (x,y) = S So(x,y), 

ZJ„<«)=/ ^..W^IAmJ «W*,*0, 



7°(i) 



I + d y b)Vip (x,y)dm(x,dy). 



Let us consider the following SDE: 
dX t = Z; /2 (X t )cfty t + Co(Xi)rft + j°(Xt)dG t , t > 0, 
X G 6, Gis continuous and increasing , 
jJ 7 °(X a )l { x. 6 e}dG. = 0,t>0, 
X = 0. 

The operators L° and T are acting on x. 

We can now state the main result of Ouknine et al. [13] (see Theorem 3.1), which is a 
generalization of a result of Tanaka [IE] (see Theorem 2.2). 



Theorem 2.1 Under assumptions (H1)-(H3), the {L £ X ,Y £ X )- diffusion process X s converges in 

law to an (L x , T x ) -diffusion X as e — > 0. 

Moreover, 

{X s , M X \G £ ) => {X, M x , G), 

where M x (resp., M ) is the martingale part of X (resp. X s ), G (resp., G £ ) is the local time 
of^{X) (resp., i){X £ ) ) at and 

" designates the weak convergence in the space C([0, T], ]R d ) x £>([0, T], x C([0, T], 
endowed with the product topology of the uniform norm for the space C([0,T],IR d ) and the 
S -topology for V([0, T],iR d ). 

It follows moreover easily from the result of Tanaka [H] that 

Lemma 2.1 Under assumptions (HI) — (H3), 

sup IE ( sup |X S £ | P + GH < oo, 

e \0<s<t J 

for any p > 1 . 

3 Reflected generalized BSDEs and weak convergence 

3.1 Preliminaries and statement of the problem 

Let us introduce some spaces: 

H 2 ' d (0,t) = 0<s<t}: IR d -valued predictable process 

such that IE J \^(s)\ 2 ds < ooj , 

S 2 (0,t) = {{ip s ; < s < t] : real-valued, progressively measurable process 



such that IE 


sup \% 


ij(s)\ 2 


< ooj 




_0<s<t 







A 2 (0,t) = {{K s ; < s < tjreal-valued adapted, continuous, , 

and increasing process such thatX(O) = 0,IE|i ; r(t)| 2 < oo} . 

In addition, we give the following assumptions: 

Let /, g : Q x ]R — >IRbea continuous function satisfying the following assumptions: 
There exist constants C > 0, p > 1, /iGlR and [3 < 0, such that for all x G 0, y, y' G IR, 
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(H.4) ^ 



' (/•«") (v - V') (/(&, y) - f(x, y')) < n\y - y'\ 2 , 
(f.n) \f(x,y)\<C(\x\P+\y\ 2 ), 

{g-i) (y - y') (g(x, y) - g(x, y')) < p\y - y'\ 2 , 



k (g.ii) \g(x,y)\ <C(\x\v+\y\ 2 ). 
Moreover, let I : — ■> IR be a continuous function and h G C 1,2 (IR + x G) for which there exist 
a constant C such that 



' (i) l(x) < C(l + \x\p), 
(ii) h(t,x) < C(l + \x\ p ) 



H.5 < 



(in) h(0,x) < l(x). 
where p is obtained in (H4). 

For each fixed (t,x) G IR + x 9, we consider (Yf, Z^, Kg) < s < t , the unique solution of 
reflected generalized BSDE 



(3.1) 



f Y s £ = l(X!) + £ f(X*, Y*)dr + f s g{Xe, Yf)dG< + Kf - K* B - £ Z^dM r x % 
Yf>h(s,Xj), 

K e is nondecreasing such that Kq = and Jq{Y £ — h(s, X £ s ))dK £ a = 0. 
Let us note that existence and uniqueness of a solution to ( 13.ip is proved in \T7\ . 

Furthermore, we also consider (Y s , Z s , K s ) < s < t the unique solution of the reflected BSDE 
f Y s = l(X t ) + Si f(X r , Y r )dr + £ g(X r , Y r )dG r - /* Z r dM r x + K t - K s 

Y s > h(s,X s ) 

K is nondecreasing such thatfQ) = and f (Y s — h(s, X s ))dK s = 0. 
Throughout the rest of this paper, we put: 



(3.2) 



Ml 



ZtdM x \ M, 



Z r dM x : < s < t 



and we consider (Y, M, K) (resp., (Y e , M £ , K 6 )) as a random element of the space D([0, t}; IR) x 
D([0,t];M) x C([0,t];M + ). D([0,t];M) denotes the Skorohod space (the space of "cadlag" 
functions), endowed with the so-called S'-topology of Jakubowski [8] and C([0, t}; IR + ) the space 
of functions of [0, t] with values in ~\R d equipped with the topology of uniform convergence. 
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3.2 Weak convergence 



The goal of this section is to prove the following result 



Theorem 3.1 Under conditions (H1)-(H4) the family (Y £ , M £ , K £ ) £>0 of processes converge 
in law to (Y, M, K) on the space D([0,t],M) x D([0,t],2R) x C([0, t\; IR+). as e -> 0. 
Moreover, 

limr £ = Y in M. 

The proof of this theorem relies on four lemmas. The first one is a simplified version of Theorem 
4.2 in (Billingsley [2], p. 25) so we omit the proof. We note also that the convergence is taken 
in law. 



Lemma 3.1 Let (U £ ) e>0 be a family of random variables defined on the same probability space. 
For each e > 0, we assume the existence of a family random variables {U £ ' n ) n<1 , such that: 

(i) For each n > 1, U £,n =^> U°' n as e goes to zero and 
(U) U°' n =^ U° asn — > +oo. 

(in) We suppose further that U £,n =^> U £ as n — > +oo, uniformly in e 
Then, U £ converges in law to U . 



For the last third lemmas, let us introduce the penalization method to construct a sequence of 
GBSDE. For each n E W, we set 

f n (x, y) = f(x, y) + n(y - h(s, x))~ 

Thanks to the result of Pardoux and Zhang [T6j, for each n G IN*, there exists a unique pair of 
^—progressively measurable processes (Y £ ' n , Z e,n ) with values in IR x IR d satisfying 



and 



where 



IE 



sup 

0<s<t 



Y £ ' n \ 2 + / \Z £ ' n \ 2 ds 



< oo 



Y r ; 



K x t ) + f fn(X £ ,Y £ ' n )dr + J' g(X £ ,Y £ ' n )dG: 
-(M £ ' n - M £ ' n ), 



M £ 



ZTdMF. 



(3.3) 



From now on, C is a generic constant that may vary from line to another. 
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Lemma 3.2 Assume (H1)-(H4). Then for any n, there exists the unique process (Y°' n M°' n ) 
such that the family of processes(Y £ ' n , M £ < n ) £>0 converges to it in D([0,T],M) x D([0,T], M), 
as e goes to 0. 



Proof. Step 1. A priori estimate uniformly in e and n 

Applying Ito's formula to the function |F/' n | 2 , we get 

|r/' n | 2 + J \Zp n \ 2 dr = \l{X £ t )\ 2 + 2 J^Y £ ' n f{X £ ,Y £ ' n )dr 

+2 f Y £ ' n g(X £ ,Y £ ' n )dG £ 



where 



-2 I Y £ ' n Z £ > n dM* £ + 2 ! Y £ ' n dK £ ' n , (3.4) 

J S J S 

K £ s ' n = n f [Y £ ' n - h(r, X £ / n )]-dr (3.5) 
Jo 

By using assumptions (H3) and (H4), the inequality J^Y £ ' n dK £ .' n < f^h(r,X £ )dK £ : n and 
Young's inequality, we get for every 7 > 

TE\Y £ ' r f + ^ \Z £ ' n \ 2 dr < (1 + CJE\X £ \ 2p ) + E ^ {(2// + 1) \Y £ > n \ 2 + C 2 \X £ \ 2q } dr 

I (2/3 + 7) \Y £ ' n \ 2 + ^ \X £ \ 2q ^ dG £ r 

h(r,X £ )dK £ ' n . 

By virtue of assumption (h.i), Young's equality and choosing 7 = \(3\ we obtain for every 8 > 
IE [\ Y s ,n \ 2 + j s \Z £ /fdr+ \(3\^ \Y £ > n \ 2 dG £ )j < C(1 + JeJ^ \Y £,n \ 2 dr) 

+7E ( sup \h(s,X £ s )\ 2 



S 



0<s<t 



+5IE (K £ ' n - K £ s ' n f . 



But 



TE(K £ t ' n - K £ ' n ) 2 < CJE (\Y £ ' n \ 2 + \l(X £ t )\ 2 + j\\K\ 2q + l^/'T)^ 

+C1E^ (\X £ T \ 2q + \Y £ ' n \ 2 )dG £ . 



f*t pt / pt 

E|Y;-"| 2 + E / \Z £ ' n \ 2 dr + E / \Y r £ ' n \ 2 dG £ r < C 1 1 + E J \Y r £ ' n \ 2 dr 



and it follows from the Gronwall's lemma that 

E ( \Y £ ' n \ 2 + [ \Z £ /fdr + / \Y £ ' n \ 2 dG £ r + \K £ t ' n \ 2 ) < C. < * <t, V n G IN*, V e > 0. 



Finally, Burkholder-Davis-Gundy inequality applying to (13 .4K gives the following: 



sup sup El sup |F/' n | 2 + / \Z £ r > n \ 2 dr+ / |Y r e ' n | 2 dG^ + \K^' n \ 2 1 < C. 
Step 2: Tightness 

We adopt the point of view of the S'-topology of Jakubovski [8]. In fact we define the 
conditional variation of the process Y £ ' n on the interval [0, t] as the quantity 



CV t (Y £ ' n ) = supE \nYtZ - Y t T/K)\\ 



where the supremum is taken over all partitions it : = t < t\ < ... < t m = t of the interval 

[o,f]. 

Clearly 

CV t (Y £ > n ) < E U\f{X £ r Xn\dr + j\g{X £ r Xn\dG £ r + \K £ t ^ , 
and it follows from step 1 and (H3) that 

sup (cV t (Y £ ' n ) + E ( sup \Y £ ' n \ + sup \M £ / l \] + \K £ ' n \ 2 ) < +oo. 

e V \0<s<t 0<s<t J J 

Let us denote 

H £ ' n = f (g(X £ } Y £ > n ) + n[Y £ ' n - h(r, X £ )}~) dG £ r . 
Jo 

By virtue of step 1 and (H3), we have 

sup ( CV t (H E ' n ) + E ( sup \H £ ' n \ j j < +oo. 

10 



Hence the sequence {(F/'™, Hp n , Mp n ); < s < t} satisfy Jakubowski tightness criterion for 
quasi-martingales under IP (see |8j), which is the same criterion as the one of Meyer-Zheng (see 

P~Q)- 

Step 3: Passage to the limit and it identification. 

After extracting a subsequence which we omit by abuse of notation, it follows by Jakubowski 
tightness criterion that (Y £,n , H £,n , M e,n ) converges in the weak sense to (Y n , H n , M n ) as s goes 
to in (-D([0, t], IR)) 3 endowed with the product of ^-topology. 

Since / is continuous, the mapping (x, y) t— > J (f(x(r),y(r)) + n[y(r) — h(r, x(r))]~) dr is con- 
tinuous from C([0, t], IR ') xD([0,t), IR), equipped respectively by topology of uniform norm and 
the S— topology. Furthermore, it follows (see Boufoussi et al.[3j) that for all < s < t 

H:= I* g(X r ,Y r °' n )dG r 

and 



o 



Y^ n = l(X t )+ f f(X r ,Y r °' n )dr+ f g(X r ,Y r °' n )dG. 

J s J s 

+n j\Y r °' n - h (r, X T ))-dr + M°' n - M°' n . 



By using a similar argument as in Pardoux [15], one can prove that M°' n and M x are Tf^ n 
martingales. Let {(Y s ,U s ),0 < s < t} denote the unique solution of the generalized BSDE 



Y" = l(X t )+J f(X r ,Y n r )dr + J g{X r X)dG r 
+n J (Y^ - h (r, X r ))~dr + - M™, 



which satisfies IE f U^d(M?)(j%)* < +oo. Set also M n s = I !T T dM* . Since F\ JT are 

Js J0 _ 

jF^-adapted, and Mf-the martingale part of X is a jF^' yn -martingale, so is M n . It follows 
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from Ito's formula for possibly discontinuous semi- martingales that 

IE IF? - Y s n \ 2 + E[M n - M n ] t - E[M n - M n ] 



= 2E / (Y r -Y r n )(f(X r ,Y r n )-f(X r ,Y r ))dr 

J s 

+21EJ (F: - Y r n )(g(X r , Y r n ) - g(X r ,Y?))dG r 
+21E J {Y n r -Y r n )(dK n -dK?) 
< (2/x + l)lE/ |F r ;-i; n | 2 dr + 2/3IE / IF™ - Y r n \ 2 dG r 



< CM I \Y n r -Y r n \ 2 dr). 

J s 



(We use the fact that (3 < 0) . 

Hence, from Gronwall lemma Y s = YJ 1 and M" = M™, < s < £. ■ 

Lemma 3.3 Assume (H1)-(H4). Then the sequence of processes (Y°' n , M 0,n , K 0,n ) n>1 con- 
verges in law to (Y, M, K) in D([0, t], iff) X D([0, t], M) x C([0, i], 2R+) endowed with the previous 
topology. 



Proof. By the same computation of step 1 and step 2 in the proof of Lemma [3721 we get 

sup IE ( sup \Y s °' n \ 2 + [ \Z^ n \ 2 dr+ [ \Y^ n \ 2 dG r + \K° t ' n f] < +oo 

and 

rO,n\ , n\/. ( js0,n\ , TTTi / Q11 _ \M0,n\2 



s 'I | r 

0<s<t ■ Jo Jo 



sup CV t (Y°' n ) + CV t (K u ' n ) + IE sup |M S U 'T < +oo. 

e V \0<s<t J J 

Further (Y n '°, K 0,n , M°' n ) satisfy Jakubowski tightness criterion for quasi-martingales under IP 
and by the similar argument as in Lemma 3.2, (Y-' n , F°' n , M°' n ), converges in law to (y, M, K\ 

as n goes to oo on D([0,t],IR 2 ) x C([0, £],IR + ) endowed with the previous topology. Finally 
by the uniqueness of the solution of the reflected generalized BSDE ( 13.21) . it follows that 

(?, M, K^j = (Y, M, K) which ends the proof. ■ 



Lemma 3.4 Assume (H1)-(H4), the family of processes (Y £ ' n , M e ' n , K e,n ) n>1 converge 
formly in e G (0, 1] in probability to the family of processes (Y e , M £ , K £ ) £e ^ x . 



um- 
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Proof. Step 1: In view of step 1 of the proof of Lemma E 



supsupIE ( sup \Y £ > n \ 2 + [ \Z £ / l \ 2 dr+ [ \Y £ ' n \ 2 dG £ . + \K £ t ' n A < +00. (3.6) 

e n \0<s<t Jo Jo J 



Step 2: There exists an jF,-adapted process {Y £ s )o< s <t such that 



lim IE / |Yl £ ' n - Yl\ 2 ds = 0. 





Notice that, for all n > 1 and (s, x, y, z) 6 [0, £] x IR d x Ex IR d , 

f n (s,x,y,z) < f n+1 (s,x,y,z). 

Therefore, by the comparison theorem for generalized BSDE |16| , we have Y £ ' n < Y £ ' n+1 , a.s. 
for all < s < t. Hence, for all < s < t 

Yr/Y £ s , a.s. 

In view of inequality (I3.6P and Fatou's lemma, 

IE ( sup \Y £ f) < C. 

\0<s<t J 

It then follows by the Lebesgue convergence theorem that 

lim IE / \Y £ > n -Yl\ 2 dr = 0. 



Step 3: lim„^ IE(sup < s < t | (Y £ ' n - h (s, X £ s ))~\ ) = a.s. 

Let |(F/' n Zf' n ), < s < t| be the unique solution of BSDE 



= l(X £ )+ / f(X £ ,Y £ ' n )dr+ / g(X £ ,Y £ ' n )dG £ r 

J s J s 

+n I (h(r,X £ )-Y £ > n )dr- f Z £ ' n dMf 

J s J s 

By using the comparison theorem, it follow for all < s < t and n > 1 that Y £,n < Y £ ' n , a.s. 
Now let v be an jF t -stopping time such that < v < t. Then, applying Ito's formula to 

ye,n e n(t-0 ) we have 

ye,n = E ^[ e -n(*-")/(xf) + j e -< s ~ v) / (X £ ,Y r £ ' n ) dr 

J v 

+ [ e- n{r ~ u) g{X £ ,Y £ ' n )dG £ + n I e' n ^h (r, X £ ) dr}. 

J V J V 
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It is easily seen that 

e~ n ^l(X £ ) +n f e- n ^h (r, X £ ) dr - l(X £ )l {l/=t} + h (t, XI) l {u<t] 

a.s. and in L 2 (fi), and the conditional expectation converges also in L 2 (Q). Moreover, 

1 



-n(r—u) 



f(X £ ,Y £ ' n )dr 



< 



t \ 1/2 



\f(X £ ,Y £ ' n )\ 2 dr 



-n(r— v) 



g(X £ ,Y £ ' n )dG £ r 



< I \g(X £ ,Y £ ' n )\ 2 dG £ 



1/2 



and 



t pt 



f(X £ ,Y r £ > n )dr + / e~ n ^g(X £ ,Yr)dG £ r 







in L 2 (Q) as n — > oo. Consequently, 

y*.« _^ /(X, £ )l {i , =t} + h {y, XI) l {v<t} in L 2 (fi) as n -> oo. 
Therefore, Fj, > h{y,X v ) a.s. 

From this and the section theorem (see Dellacherie et Meyer [4] page 220), the previous in- 
equality hold for all s e [0, t] a.s. Then 

(Y s n - h{s,X s ))~ \ 0, < s < t, a.s., 

and from Dini's theorem the convergence is uniform in t. Since (Y s n — h(s,X s ))~ < (h(s,X s ) — 
F,°) + < \h(s,X a )\ + \Y S \, the result follows by dominated convergence. 

Step 4: lim n ^ IE (sup 0<s<t \Y £ ' n — Y s \ 2 ) = and there exists jF r adapted processes (Z s ) < s <t 
and (K s ) < s < t such that 



limE 

n— > 



Z\ 2 + sup \K £ s ' n -K s \ 2 | = 0. 

0<s<t 



Indeed, by Ito's formula and taking expectation in both sides, we have for all n > m > 1 
and < s < t, 



|y/' n - Y s £ > m \ 2 + I \Z £ r > n - Z £ r > m \ 2 dr + 2\{3\ I \Y £ > n -Y £ > m \ 2 dG 



< 2fj, J \Y £ ' n - Y £ ' m \ 2 dr + 2 J (Y £ ' n - h (t, X £ )YdK\ 



r 



■ 2 / ( ) - // (r, X £ r )YdK £ : m - / (Y £ > n - Y £ > m ) (Z £ .' n - Z £ : m )dMf 
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(3.7) 



Therefore 



IE [ \Z £ r ' n - Z E r ' m \ 2 dr < CJE I \Y r E ' n - Y r £ ' m \ 2 dr + 2JE f {Y r £ ' n -h{r,X £ r ))-dK £ r 

J s J s J s 

+2E I {Y £ > m -h{r,X £ .))-dK £ ' m . 

J s 

Furthermore, Holder's inequality provide 



IE / (Y £ > n - h (r,X £ ))-dK £ ' n < 



IE sup \(Y?> n -h(s,X,))- 



and 



IE J (Y £ ' p - h (r, X £ ))~dK £ ' n < 



0<s<t 



IE sup \(Y s £ ' m -h(s,X s )) 

0<s<t 



q 1/2 



[ E | X m|2]l/2 (3^ 



-|2 



1/2 



[E|K™| 2 ] 1/2 . (3.9) 



By using the result of Step 1 and Step 2, the right hands sides of the above inequalities (13.81) 
and (13.91) go to zero asm,n-> oo. It follows that 



7S,n ye,m\ 



dr -> 



o 



as n, m — > oo, and there exists a process Z e H ' (0,t) such that 

nt 

IE / \Z £ ' n -T r \ 2 dr -> 0, asrwoo. 
Then, by the Burkholder-Davis-Gundy inequality, we get 

IE( sup \Y £ ' n - Y £ ' rn \ 2 < 4IE / (Y n - S s ydK™ + 4IE [ (Y m - S s )~dK r s 

0<s<t Jo JO 

+CIE f \Z £ ' n - Z £ ' m \ds. 
Jo 

It follows that IE(sup 0<s<t \Y £,n — Y £ ' m \ 2 — > asn, m — > oo. Therefore, 



We deduce that 



ie<; sup |F/'"-r^| 2 

0<s<t 



IE ( sup |if s e ' n - K £ s ' m \ 

0<s<t 



as 72 — > oo . 



asn, m — > oo. 
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Hence there exists a process (-^g)o<s<t £ ^ 2 (0,£) such that 

E( sup |K S E '"-^| 2 ] ->0asra->oo. 



0<s<t 



Step 5: The limiting process {Y £ , Z e , if 6 ) is the solution of the reflected generalized BSDE 
Obviously the process (Y s , Z s , K s )o< s <t satisfies 



y £ s = i+ I f(r,X s r ,Y £ r ,Z £ r )ds+ I g(r,X £ r ,Y £ r )dG e r + K £ t -K £ s - [ Z £ r dM?\ 0<s<t. 

J s J s J s 



Since (Y e ' n , K £,n ) tends to (Y £ , K £s j in probability uniformly in t , the measure dK £,n converges 

to dK £ weakly in probability, so that f (Yf' n — h (r, X £ ))dKp n — > f Q (Y e r — h (r, X £ ))dK £ r as 

n — > oo. Obviously, fn(Y e r — h (r, X £ ))dK £ r > 0, while, on the other hand, for all n > 0, 

Jo(y r e ' n - /i (r, X £ r ))dK £ r > n < 0. 
Hence, 



/ (F^-/i(s,X s e ))d^ = 0. a.s. 



Consequently, (Y £ s , Z 5 ,, K £ s )o< s < t is the solution of the reflected generalized BSDE (£, /, g, /i) and 
uniqueness of this equation allows that Y = Y £ , Z = Z £ and K = K £ which ends the proof. 



Proof of Theorem [37T1 

Put Theorem 2.1, Lemma [321 Lemma l3~3l Lemma [3T41 and the Proposition 13.11 together prove 
that {Y £ , M e , K £ ) converges in law to (Y, M, K) in D([0,T],IR) x L>([0, T], IR) x C([0, T], M + ), 
where the first factor is equipped with the S'-topology, and the two last factors with the topology 
of uniform convergence. In particular 

Y £ = IE il (X £ ) + jf / {XI Y £ ) da + J*g {X £ s , Y £ ) dG £ s + K £ t 

converge to 

r = iE + ^ f(x a ,Y s )ds + jf ^(x S) y s )dG s + ^| 

in IR. Indeed, since YJf is deterministic, Yq = TE(B e ) where 

B £ = I (X £ ) + f f {XI Y £ ) ds+ fg {XI Y £ ) dG £ s + K £ t . 
Jo Jo 
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In view of assumptions (H1)-(H4), we have, 

supEl^l 2 < C. 



£ 



Since B £ converge in law as e goes to 0, toward 



B = l(X t )+ [ f(X s ,Y s )ds + [ giX.MdG. + Kt 



ti 



Jo Jo 
the uniformly integrability of B £ implies that 



limlE(S e ) = IE (hmB £ ) . 



This mean that Y £ converges to Y = JE(B). 
■ 

4 Homogenization of the reflected viscosity solutions of 
scalar PDE with non linear Neumann boundary condi- 



The result of the above section permits to us to deduce weak convergence of a sequence (Y e ) 
of the solutions of reflected generalized BSDEs from weak convergence of the sequence (X £ ). 
In a sense, we deduce, from the probabilistic proof of the convergence of linear PDE with 
linear boundary condition due to Tanaka [18], a probabilistic proof of convergence of reflected 
semi-linear PDEs with non linear Neumann boundary condition. 

For each x e 0, let {{X £ s ' x , G £ s ' x ), s > 0} denote the solution of the SDE (IQjl . Let I, f, g, h 
be as in the previous section. For each (t,x) G IR + x 0, let 



where Y £ denote the solution of the reflected generalized BSDE (13.11) considered in the previous 
section (which of course depends on the starting point x of X £ , and the final time t for the 
reflected generalized BSDE). The continuity of u £ follows by the standard arguments (see [T7] ) 
from the continuity of the mapping x i— ► X £,x in probability. One can show as in Ren and Xia 
[17] that u £ is the viscosity solution of the semi-linear parabolic PDEs 



tion 



u £ (t,x) = Y 



r£,X 
5 



min [u £ (t, x) — h (t, x) , 



9m e 

at 



(t, x) - U x u £ (t, x)-f (x, u £ (t, x))) = 0, i 



f (t,x) e IR+ x 0, 



< 



(4.1) 



Y £ u £ (t, x) + g (x, u £ (t, x)) — 0, i 



f (t, x) e IR + x dQ, 



u £ (0, x) = I (x) , x E 0, 
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where U x and T e x are the operators defined by (12. ip and (12.21) respectively. 
Let u be the solution of the homogenized PDE 
min {u (t, x) — h (t, x) , 

-f (t, x) - L° x u (t, x)-f (x, u (t, x))) = 0, if (t, x) G H + x 6, 

(4.2) 

Tu (t, x)+g (x, u (t, x)) = 0, if (t, x) G IR + x d&, 
u (0, x) = I (x) , x e 0, 

where L° x and T° are defined in ( 12.41) . Recall again the previous argument, u(t, x) = Yq is 
the viscosity solution of ( 14. 2ft . where F 21 denote the solution of the reflected generalized BSDE 
(13.2j) . which depend on the starting point of JT X . 

Theorem 4.1 Under the assumptions (H1)-(H4), for each (t,x) G iR + x 9, T/zen u £ (t,x) 
converges to u (t, x) as e goes to 0. 

Proof. Since u £ (t,x) = Yq' x and u(t, x) = Yq, the result follows from Theorem 3.1 ■ 
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